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Abstract 

The mass and kinetic energy distribution of nuclear fragments from thermal neutron induced 
fission of ^35^ have been studied using a Monte-Carlo simulation. Besides reproducing the 
pronounced broadening on the standard deviation of the final fragment kinetic energy distribution 
(cje(m)) around the mass number m = 109, our simulation also produces a second broadening 
around m = 125, that is in agreement with the experimental data obtained by Belhafaf et al. 
These results are consequence of the characteristics of the neutron emission, the variation in the 
primary fragment mean kinetic energy and the yield as a function of the mass. 

KeywordsrMonte-Carlo; neutron induced fission; '^^^U; standard deviation. 

Mediante la simulacion con el metodo Monte-Carlo fue estudiado la distribucion de masas y 
energia cinetica de los fragmentos de la fision inducida por neutrones termicos del ^^^[/. Ademas 
de reproducir el ensanchamiento pronunciado en la desviacion estandar de la distribucion de la 
energia cinetica de los fragmentos finales {ae{m)) alrededor del numero masico m = 109, nuestra 
simulacion tambien produce un segundo ensanchamiento alrededor de m = 125, en concordancia 
con los datos experimentales obtenidos por Belhafaf et al. Estos resultados son consecuencia 
de las caracteristicas de la emision de neutrones, la variacion de la energia cinetica media y el 
rendimiento de los fragmentos primarios en funcion de la masa. 

Descriptores: Monte-Carlo; fision inducida por neutrones; ^^^C/; desviacion estandar. 

PACS: 21.10.Gv; 25.85.Ec; 24.10.Lx 
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I. INTRODUCTION 



Since the discovery of the neutron-induced fission of uranium by Hahn and Strassmann 
in 1938 jl|, much effort has been made to understand the processes involved in it and to 
measure the relevant fission parameters. Nowadays several aspects of heavy nuclei fission 
seem to be clarified. Meitner and Frisch suffffested a theoretical explanation based on a nu- 
clear liquid-drop model [2j , and, over the past 30 years the model has provided considerable 
insight into nuclear structure [3]. It is known that the de-excitation by fission of heavy nuclei 
depends of the quantum properties of the saddle point and of the associated fission barrier. 
The detection of fission isomers has been interpreted by the secondary well in the fission 
barrier ^. The nascent fragments begin to be formed at the saddle point, then the system 
falls down to the fission valley (energetically preferred paths to fission) and ends at the scis- 
sion configuration where fragments interact only by Coulomb force. Moreover, at scission, 
the fragments have acquired a pre-scission kinetic energy. Over the fission valley, the system 
could be described by collective variables (such as deformation, vibration, rotation, etc.) 
and intrinsic variables (such as quasi-particles excitations). Nevertheless, the dynamics of 
the fission processes are not yet completely understood . In particular, it is neither known 
the nature of the coupling between the collective and intrinsic degrees of freedom during the 
descend from the saddle to scission, nor known how it does arise. The physics problem of 
the description of the fission fragment mass and kinetic energy distributions is very closely 
related to the topological features in the mult i- dimensional potential energy surface . In 
the low-energy fission, several final fragment characteristics can be explained in terms of 
a static scission model of two coaxial juxtaposed deformed spheroidal fragments, provided 
shell effects, affecting the deformation energy of the fragments. These shell effects correc- 
tions, determined by the Strutinsky prescription and discussed by Dickmann et al. [3] and 
Wilkins ^, subsequently generate secondary minima in the total potential energy surface 
corresponding to fragments having some particular neutron or proton shell configurations. 
If the final fragment characteristics were governed by the properties of the fragments them- 
selves, a basic argument in any statistical theory, one would then expect an increase in 
the width of the kinetic energy distribution curve for fragment masses A, having the above 
mentioned special neutron or proton shell arrangements. In order to address this question, 
the fission parameters of the primary fragments (pre-neutron emission) have been the most 
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FIG. 1: Thermal neutron induced fission of ^^^C/. Standard deviation of the final fragment kinetic 
energy distribution as a function of the final mass m, as a result of Monte-Carlo simulation (A), 
and experimental data (•). Both from Q]. 

studied are the mass yield (y(y4)) and the kinetic energy [Ei^A]) distribution. Nevertheless, 
direct measurements can only be carried out on the final fragments (post neutron emission) 
mass yield Yim) and kinetic energy (e(m)). Therefore it is crucial to find out what is the 
relation between the primary and the final kinetic energy distributions, as well as the re- 
lation between the Yi^A) and Yim) curves. For thermal neutron induced fission of ^^^f/, 
which in fact is the fission of excited ^ssfj ^236^*^ formed by neutron absorption by ^^^t/, 
the e(m) distribution was experimentally determined by Brissot et al. 9|]. This distribution 
was represented by the mean value of kinetic energy e and the standard deviation (SD) of 
the kinetic energy o"e as function of the final mass m. As seen in Fig. [T]the plot of both the 
measured values and the results of a Monte-Carlo simulation of CTg from a primary distribu- 
tion E{A) without broadenings, shows one pronounced broadening around m ^ 109. This 
Monte-Carlo simulation result suggests that the broadening does not exist on the primary 
fragment kinetic energy as a function of the primary fragment mass. In a latter experiment, 
Belhafaf et al. [l^, repeated the experiment of Brissot et al. for neutron induced fis- 
sion of ^^^[/, obtaining a second broadening around m ^ 125 (see Fig. 2). A Monte-Carlo 



4 



> 



Q 
CO 




110 120 

Fragment Mass 

FIG. 2: Thermal neutron induced fission of '^^^U. Simulated standard deviation of the final frag- 



ment kinetic energy distribution as a function of the final mass m (A), from Ref. 

□ 

reproduce the experimental broadening around m = 125 taken from Ref. [lOl ]. 



does not 



simulation made by these authors, from a primary distribution of E{A) without a broade- 
ning, reproduced the experimental broadening on CTg at m = 109, but failed to reproduce 
the broadening around m = 125. They suggested that this broadening must exist in the 
primary fragment kinetic energy {E{A)) distribution, and accordingly they fitted their ex- 
perimental data from a distribution with a broadening around ^4=126. In this paper, we 
present new Monte-Carlo simulation results for thermal neutron induced fission of ^^^[/. We 
compute both the mass and kinetic energy of the primary and final fission fragments, and we 
show that the broadenings on the cXe curve around the final fragment masses m = 109 and 
m = 125 can be reproduced without assuming an adhoc initial structure on aE{A) curve. 

II. MONTE CARLO SIMULATION MODEL 

A. Fragment kinetic energy and neutron multiplicity 



In the process of thermal neutron induced fission of ^^^[/, the excited composed nucleus 
236 fj* jg formed first. Then, this nucleus splits in two complementary fragments having Ai 
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FIG. 3: Thermal neutron induced fission of '^^^U . Simulation results for the primary (A) and final 
(O) mass yields are presented together with experimental data (•), taken from Ref. ll| 



and A2 as mass numbers, and E\ and E2 as kinetic energies, respectively. 

Using relations based on momentum and energy conservation, the total kinetic energy of 
complementary fragments results 

TKE = Ei + E, = ^'I^' Ei. (1) 

The total excitation energy is given by 

TXE = g - e„ - TKE, (2) 

where Q is the difference between fissioning nucleus mass and the sum of two complementary 
fragments masses, and e„ is the separation neutron energy of "^^^U. Using equation ([1]) in 
([2]) and taking into account that Ai + A2 = 236 results 

236 

TXE = Q + en tE, (3) 

^ " 236 -yl ' ^ ^ 

where A and E are the mass number and kinetic energy, respectively, of one of the two 
complementary fragments. It is reasonable to assume that the excitation energy of one 



6 



complementary fragment (i?*) is proportional to the total excitation energy, then, 

236 

E* oc TXE = Q + e„ - ^^^^rjE, (4) 

and that the number (z/) of neutrons emitted by a fragment is proportional to its excitation 
energy, i.e. 

u oc E*. (5) 
From relations (jlj) and ([5]) one derives a linear relation between u and E: 

u = a + bE. (6) 

Taking into account that there is no neutron emission z/ = for fragments having the 
maximal kinetic energy (Emax) and assuming that for the average value of fragment kinetic 
energy ly = u, the relation results 

E - E 

Let be the parameter /3 define the maximal value of kinetic energy by the relation 

Emax ^ ^ ~^ 

Then, the relation ([71) may be expressed as 

u = u{l-P{^^))- (9) 
Because the neutron number is integer, it will be defined as the integer part of ([9]), i.e. 

Jt; 

N = Integer part of(a + - p{ ))), (10) 

where a is used to compensate the effect of the change from a real number u to an integer 
number A^. 



B. Simulation process 

In our Monte Carlo simulation the input quantities are the primary fragment yield (Y), 
the average kinetic energy (E), the standard deviation of the kinetic energy distribution 
{(Je) and the average number of emitted neutron (u) as a function of primary fragment mass 
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(A). The output of the simulation for the final fragment are the yield (Y), the standard 
deviation of the kinetic energy distribution (cr^;) and the average number of emitted neutron 
(z/) as a function of final fragment mass m. 



For the first simulation, we take Y from Ref. 
et al. 



121 1. and E from Ref. 



Ill ], p from experimental results by Nishio 



101 ]. The first standard deviation ge curve is taken without 
any broadening as function of A. Then, we adjust Y{A), E{A) and aE{A) in order to 

get y(m), z/, e(m), (Je(m) in agreement to experimental data. 

In the simulation, for each primary mass A, the kinetic energy of the fission fragments is 
chosen randomly from a Gaussian distribution 



P{E) = — — exp 
\J2i\oe 



{E-Ef 



(11) 



where P{E) is the probability density of energy with mean value E and standard deviation 
(Te- 

For each E value, the simulated number of neutrons N is calculated with the relation 
f fTOj) . The final mass of the fragment will be, 

m = A-N. (12) 

Furthermore, assuming that the fragments loose energy only by neutron evaporation and 
not by gamma emission or any other process, and neglecting the recoil effect due to neutron 
emission, then the kinetic energy e(m) of the final fragment will be given by 

e(m) = (1 - ^)E. (13) 

With the assemble of values corresponding to m, e and A^, we calculate Y{m), e(m), cre(m) 
and z/(m). 

On the other hand, to obtain an acceptable statistics during the simulation, we have 
considered a total number of fission events of ^ssf/ ^-j-^g order of 10*^. At the same time, 
we have used the Box-MuUer method to generate the random numbers with the required 
normal distribution [l^, and have computed the SD of all the relevant quantities by means 
of the following expression which for e(m), read as 

o {rn)= \ \ e (m), 14) 

where e(m) is the mean value of the kinetic energy of final fragments with a given mass m, 
and Nj{m) is the number of fission events corresponding to that mass. 
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FIG. 4: Thermal neutron induced fission of ^^^C/. Mean kinetic energy of the final fragment (0) 
and the mean kinetic energy of the primary fragments A, as a result of simulation in this work, to 
be compared to experimental data (•) taken from Ref. 
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III. RESULTS AND DISCUSSION 



The simulated final mass yield curve Y{m) and the primary mass yield curve Y{A) are 
illustrated in Fig. [31 As expected, due to neutron emission, the Y{m) curve is shifted from 
Y{A) towards smaller fragment masses. As stated in sect. [Tll the primary kinetic energy 
{E{A)) is generated from a Gaussian distribution, while the final kinetic energy (e(m)) is 
calculated through Eq. (1131) . The plots of the simulated mean kinetic energy for the primary 
and final fragments as function of their corresponding masses, are shown in Fig. H] . In 
general, the simulated average final kinetic energy curve as a function of final mass (e(m)) 
have roughly a shift similar to that of Y{m) curve, and a diminishing given by relation 
f|T3|) with N = v. The exceptions of this rule are produced in mass regions corresponding 
to variations of the slope of Y{A) or E{A) curves, for example for A = 109, A = 125 
and A = 130. Furthermore, Fig. \5\ displays the standard deviation of the kinetic energy 
distribution of the primary fragments and the standard deviation of the kinetic energy of the 
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FIG. 5: Thermal neutron induced fission of U. Standard deviation of final fragment kinetic 
energy distribution (0) and standard deviation of primary fragments kinetic energ y d istribution 
(A), as simulated in this work, to be compared to experimental data (•) from Ref. id ]. 

final fragments (cre(m)). The plots of c'"e(m) reveal the presence of a pronounced broadening 
around m = 109, and a second broadening is found around m = 125, in a mass region 
where there are variations of the slopes of Y{A) or E{A) curves. There is no experimental 
data around m = 130. Nevertheless, if one takes the experimental value cTp = 3.9MeV for 

n 

m = 129 from Ref. M and one puts it on Fig. \5\ the beginning of another broadening for 
m = 130 is suggested. 

These results were obtained with a simulated primary fragment kinetic energy distribution 
(see Fig. OA) without broadenings in the range of fragment masses A from 90 to 145. If 
one simulates an additional source of energy dispersion in as, without any broadening, no 
broadening will be observed on ae- 

Both the shape and height of the broadenings of (Jeifn) are sensitive to the value of pa- 
rameter a and 13 appearing in Eq. ( fTOi) . A higher value of a will produce a larger broadening 
of SD. The effect of (3 on broadening depends much on mass region. For the region m = 109, 
a higher value of /9 will produce a larger broadening of SD. The simulated results for cre(m) 
presented in Fig. [5] were obtained with a = 0.62 and /?=0.35. 
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FIG. 6: The average number of emitted neutrons from fission of ^35^/ . function of the primary 
fragment mass A (A), as a function of final fragment mass (0) both as result of simulation and 
experimental data (•), taken from Ref. [2] 

The simulated average number of emitted neutron z^(m) curve is shifted from iy{A) in a 
similar way as y(m) relative to y(v4)(see Fig. [6]). 

The presence of broadenings about m = 109 could be associated with neutron emis- 
sion characteristics (approximately u = 2) and a very sharp fall in kinetic energy from 
E =100 MeV to E =85.5 MeV, corresponding to ^=109 and ^=111, respectively. The 
second broadening is produced by a discontinuity of the curve E{A) between A =126 to 
A =125, which is necessary to reproduce a similar discontinuity between m =125 to m 
=124. We give emphasis to the shape of Ue which increase from m = 121 to m = 125 and 
it decreases from m = 125 to m = 129 as occurs with experimental data. 



IV. CONCLUSION 



Using a simple model for the neutron emission by fragments, we have carried out a 
Monte-Carlo simulation for the mass and kinetic energy distributions of final fragments from 
thermal neutron induced fission of "^^^U. In comparison with the primary fragments, the final 



11 



fission fragments liave eroded kinetic energy and mass values, as mucli as to give rise to tlie 
appearance of broadenings in tlie standard deviation of the final fragments kinetic energy as 
a function of mass cre(m) around m = 109 and m = 125 respectively. These broadenings are 
consequence of neutron emission and variations on slopes of primary fragments yield {Y{A)) 
and mean kinetic energy E{A) curves. From our simulation results, another broadening, 
around m = 130, may be predicted. 
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